#tesis

import numpy
#import 

#Entradas
#G (alfa,beta)	: (W/m2) Irradiancia global sobre el plano del generador

Gdm_avg=(6.35,6.49,6.30,5.76,4.70,3.72,3.61,4.05,4.98,5.59,6.10,6.35)

#Isc_std
#Ta	:temperatura ambiente 
#TONC
#Voc_std
#FFstd
#Nmp
#Ncp
#Nms
#Ncs

#Calculadas
#Ics :Corriente de corto circuito 
#Tc
#Voc
#Vt
#voc
#FFo
#rs
#a
#b
#Vmax
#Imax
#Pmax
#Vmaxg
#Imaxg
#Pmaxg
#-----------------

#Calculo de irradiancia en el plano  del generador


#Radiacion directa y difusa  horizontales diarias medias mensuales 

#Estimacion de radiacion directa o difusa en intervalos de una hora


#Estimacion de radiacion Global (Gb+Gd+Gr) en intervalos de una hora en una superficie inclinada


#---------------
#Corriente de corto circuito
Isc=G*Isc_std/1000

#Temperatura de celula
Tc=Ta+(TONC-20)/800*G

#Tension circuito abierto de la celula
Voc=Voc_std-0.0023*(Tc-25)

#Voltaje termico
Vt=0,025*(Tc+273)/300

#tension de celula normalizada
voc=Voc/Vt

#Factor de forma de una celula ideal sin considerar la resistencia serie
FFo=(voc-numpy.log(voc+0.72))/(voc-1) 

#Resistencia normalizada
rs=1-FFstd/FFo

#Tension y corriente de celula en el punto de maxima potencia
a=voc+1-2*voc*rs
b=a/(1+a)

Vmax=Voc*(1-b/voc*numpy.log(a)-rs*(1-a**-b))
Imax=Isc*(1-a**-b)

#Maxima Potencia de celula
Pmax=Vmax*Imax

#Para el generador
Vmaxg=Vmax*Nms*Ncs
Imaxg=Imax*Nmp*Ncp
Pmaxg=Pmax*Nmp*Ncp*Nms*Ncs


